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I. Introduction 

The search for Higgs bosons is a crucial test of Supersymmetry (SUSY) which can be 
performed with the present and the next generation of accelerators. A precise prediction 
for the production and decay processes of the Higgs bosons in terms of the relevant SUSY 
parameters is necessary in order to determine the discovery and exclusion potential of the 
Tevatron [1-3], and for physics at the LHC [4-7] and the ILC [8-12]. A precise predic- 
tion of Higgs-boson properties is also required as input for (current) electroweak precision 
analyses, see e.g. Refs. [13,14] 

Since the soft SUSY-breaking parameters can in general be complex, it is necessary to 
take into account the effects of complex phases. In the Minimal Supersymmetric Standard 
Model with complex parameters (cMSSM) Higgs physics is affected by complex param- 
eters entering via loop corrections. In particular, these are the Higgs mixing parameter, 

II, the trilinear couplings, Af, f = t.h.r, . . ., and the gaugino mass parameters Mi, M2, 
M3, where jMsj = rrig (the gluino mass). As a consequence, once loop corrections are 
taken into account the neutral Higgs bosons are no longer CP-eigenstates, but mix with 
each other [15], 



As tree-level input parameters in the Higgs sector (besides the gauge couplings and the 
Z boson mass) it is convenient to choose the mass of the charged Higgs boson, Mfj± , and 



{h,H,A) hi,h2,hs ■ 



(1) 
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the ratio of the two vacuum expectation values, tan /3. 

In order to obtain the prediction for the Higgs masses beyond lowest order, the poles of 
the Higgs propagators have to be determined. Since the propagator poles are located in the 
complex plane, we define the physical mass of each particle according to the real part of 
the complex pole. 

Neglecting the mixing with the Goldstone bosons and the (longitudinal parts) of the 
Z boson (see Ref. [16] for details) one can write the propagator matrix of the neutral Higgs 
bosons /i, il, A as a 3 X 3 matrix, A/ijja(p^)- It is related to the 3 x 3 matrix of the 
irreducible vertex functions by 

i^hHA{p') = -{thHA(p^'))\ (2) 

where 

r,»HA(p') = i[p'l-M„(p2)] , (3) 

Mn{p')^\ -±hHip^) ml-tHH{p^) -tHA{p^) |- (4) 

-thAip") -^ha{p^) ml-±AAip 



Inversion of t'hHAip^) yields for the diagonal Higgs propagators (i = h, H, A) 

^"^P') = zlveff/ 2^ ' 

p^ -m| +I]^f 

where Ahhip"^), Ahh(p^), Aaa(p^) are the (11), (22), (33) elements of the 3x3 matrix 
^hHA{p^), respectively. The structure of eq. (5) is formally the same as for the case 
without mixing, but the usual self-energy is replaced by the effective quantity £|f (p^) 
which contains mixing contributions of the three Higgs bosons. It reads (no summation 
over i,j,k) 



Sf(p^) = S,,(p^)-i- 



r,,(p2)r,,(p2)_r2^(p2) 

(6) 



where the f ^ (p^) are the elements of the 3 x 3 matrix f ^//^(p^) as specified in eq. (3). 

For completeness, we also state the expression for the off-diagonal Higgs propagators. 
It reads (i ^ j, no summation over k) 

Ay (p2 ) = . ^'■'I'^'''' ^^^'1^^' . - ^ , (7) 

where we have dropped the argument p^ of the Tij {p^) appearing on right-hand side for 
ease of notation. 

The complex pole A^^ of each propagator is determined as the solution of 

Ml-ml + m{M^i) = Q- (8) 
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Writing the complex pole as 



M"^ = M^ - iMT, 



(9) 



where M is the mass of the particle and T its width, and expanding up to first order in T 
around yields the following equation for Mf, 



The short-hand notation f'{p'^) = d f (p^) / {d p^) has been used, and M, denotes the loop- 
corrected mass, while mj is the lowest-order mass (i = h, H, A). 

While the Higgs-boson masses Mf can in principle directly be determined from eq. (10) 
by means of an iterative procedure (since appears as argument of the self-energies in 
eq. (10)), it is often more convenient to determine the mass eigenvalues from a diagonal- 
ization of the mass matrix in eq. (4). In the code FeynHiggs 2 . 5we perform a numerical 
diagonalization of eq. (4) using an iterative Jacobi-type algorithm [17J. The mass eigen- 
values Mi are then determined as the zeros of the function nf{p^) — p^, where Hiip"^) is 
the ith eigenvalue of the mass matrix in eq. (4) evaluated at p"^. Insertion of the resulting 
eigenvalues Mi into eq. (10) verifies (to C(r)) that each eigenvalue indeed corresponds to 
the appropriate (complex pole) solution of the propagator. We define the loop-corrected 
mass eigenvalues according to 



In our determination of the Higgs-boson masses we take into account all imaginary 
parts of the Higgs-boson self-energies (besides the term with imaginary parts appearing 
explicitly in eq. (10), there are also products of imaginary parts in ReS^?(M?)). The ef- 
fects of the imaginary parts of the Higgs-boson self-energies on Higgs phenomenology 
can be especially relevant if the masses are close to each other. This has been analyzed 
in Ref. [18] taking into account the mixing between the two heavy neutral Higgs bosons, 
where the complex mass matrix has been diagonalized with a complex mixing angle, re- 
sulting in a non-unitary mixing matrix. The effects of imaginary parts of the Higgs-boson 
self-energies on physical processes with s-channel resonating Higgs bosons are discussed 
in Refs. [18-20]. In Ref. [18] only the one-loop corrections from the t/i sector have been 
taken into account for the H-A mixing, analyzing the effects on resonant Higgs production 
at a photon collider. In Ref. [19] the full one-loop imaginary parts of the self-energies have 
been evaluated for the mixing of the three neutral MSSM Higgs bosons. The effects have 
been analyzed for resonant Higgs production at the LHC, the ILC and a photon colhder 
(however, the corresponding effects on the Higgs-boson masses have been neglected). In 
Ref. [20] the t/b one-loop contributions (neglecting the t/h corrections) on the H-A mix- 
ing for resonant Higgs production at a muon collider have been discussed. Our calculation 
incorporates for the first time the complete effects arising from the imaginary parts of the 
one-loop self-energies in the neutral Higgs-boson propagator matrix, including their effects 
on the Higgs masses and the Higgs couplings in a consistent way. 

As described above, the solution for the Higgs-boson masses in the general case where 
the full momentum dependence and all imaginary parts of the Higgs-boson self-energies 




(10) 



Mh^ < Mh, < Mh^. 



(11) 
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are taken into account is numerically quite involved. It is therefore of interest to consider 
also approximate methods for determining the Higgs-boson masses. Instead of keeping the 
fuU momentum dependence in eq. (4), the "p^ on-shell" approximation consists of setting 
the arguments of the self-energies appearing in eq. (4) to the tree-level masses according 

toiij = h, H, A) 

on-shell approximation: Eii(p^) ^ Sii(m^) (12) 
Sy(p2)^S,,((m,2+m2)/2). 

In this way the Higgs-boson masses can simply be obtained as the eigenvalues of the 
(momentum-inpendent) matrix of eq. (4). The "p^ on-sheU" approximation has the benefit 
that it removes all residual dependencies on the field renormaUzation constants that cannot 
be avoided in the iterative procedure. 

Instead of setting the momentum argument of the renormalized self-energies to the 
tree-level masses, in the "p^ = 0" approximation the momentum dependence of the self- 
energies is neglected completely (i, j = h, H, A), 

p^ = approximation: Sii(p^) Sii(O) (13) 
Ei,(p2)-.S,,.(0). 

In the "p^ = 0" approximation the masses are identified with the eigenvalues of Mn (0) 
(see eq. (4)) instead of the true pole masses. This approximation is mainly useful for 
comparisons with effective-potential calculations. The matrix Mn(0) is hermitian by con- 
struction. 



2. Amplitudes with external Higgs Bosons 

In evaluating processes with external (on-shell) Higgs bosons beyond lowest order one has 
to account for the mixing between the Higgs bosons in order to ensure that the outgoing 
particle has the correct on-shell properties such that the S matrix is properly normalized. 
This gives rise to finite wave-function normalization factors.^ For the case of 2 x 2 mixing 
appearing in the MSSM with real parameters (rMSSM) for the mixing between the two 
neutral CP-even Higgs bosons h and H, which is analogous to the mixing of the photon 
and Z boson in the Standard Model, the relevant wave function normalization factors are 
well-known, see e.g. Refs. [23,24]. An amplitude with an external Higgs boson, i, receives 
the corrections {i,j = h, H, no summation over 

(Vi + ZijVj) {i ^ j) , (14) 



^The introduction of these factors can in principle be avoided by using a renormalization scheme 
where all involved particles obey on-shell conditions from the start, but it is often more convenient 
to work in a different scheme like the DR scheme for the field renormalizations [22]. 
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where the Ti^ denote the one-particle irreducible Higgs vertices, and 



Zi = 



l + ReE^(p2)-Re 



Zii — — 



(15) 



(16) 



As before rrij denotes the tree-level mass, while is the loop-corrected mass. 

In the case of the cMSSM, the formulas above need to be extended to the case of 3 x 3 
mixing. A vertex with an external Higgs boson, i, has the form (with k all different, 
i,j,k = h, H, A, and no sunnmation over indices) 



(17) 



where the ellipsis represents contributions from the mixing with the Goldstone boson and 
the Z boson, see Ref. [16] for more details. The finite Z factors are given by 



1 



(ReSf)'(M2)' 



(18) 



1 + 



Zi-i — 



p^=Mf 



Mi - ml + ±kk{Mi)) - E,fc(M2)S,,(M2) 



%(M2) - - m| + S,,(M2)) -ml + t,,{Mi)) ' 



(19) 



where the propagators Aji(p^), Ajj(p^) have been given in eqs. (5) and (7), respectively. 
Using eq. (17) with Z factors specified in eqs. (18), (19) and adding to this expression 
the mixing contributions of the Higgs bosons with the Goldstone bosons and the gauge 
bosons (see Ref. [16] for more details) yields the correct normalization of the outgoing 
Higgs bosons in the S matrix. 

For convenience we define a matrix Zn based on the wave function normalization fac- 
tors. Its elements are given by (with Zu = 1, i,j = h,H,A, and no summation over 
i) 



(Z„)i 



Zi Zij . 



(20) 



Performing a re-ordering of the lines of Zn such that they correspond to the mass ordering 
of eq. (11) results in the matrix Z„. A vertex with an external (on-shell) Higgs boson hi is 
then given by 



(21) 



where the ellipsis again represents contributions from the mixing with the Goldstone boson 
and the Z boson. 
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3. Effective couplings 



In a general amplitude with internal Higgs bosons, the structure describing the Higgs part 
is given by Tj, where the Tij are as above the one-particle irreducible Higgs 

vertices, and the propagators are given in eqs. (5) and (7). For phenomenological 
analyses it is often convenient to use approximations of improved-Bom type with effective 
couplings incorporating leading higher-order effects. There is no unique prescription how 
to define such effective coupling terms. One possibility would be to consider the matrix 
Zn, defined through eqs. (20)-(21), as mixing matrix. The elements of the matrix Zn, 
however, are in general complex, so that the Zn is a non-unitary matrix. Therefore it 
cannot be interpreted as a rotation matrix. If one wants to introduce effective couplings by 
means of a (unitary) rotation matrix, it is necessary to make further approximations. 

A possible choice leading to a unitary rotation matrix is the "p^ = 0" approximation, 
which is used in the effective potential approach. As before, we first consider the case of 
2x2 mixing relevant for the rMSSM. In the "p^ = 0" approximation defined in eq. (13) 
the momentum dependence in the renormalized self-energies is set to the respective tree- 
level Higgs boson masses, so that the derivative in eq. (15) acts only on the term in the 
propagator factor. In this Umit Zi simplifies to [25,26] 

= approximation, 2x2 mixing: Zj = . (22) 

For the mixing between the neutral CP-even Higgs bosons /i, U this yields Zh = Zh = 

cos^ Aa. This corresponds to an effective coupling approximation where the tree-level 
mixing angle a appearing in the couplings of the neutral CP-even Higgs bosons is replaced 
by Q!eflF = a + Aa [25,26]. 

It is easy to verify that for the 3x3 mixing case eq. (18) in the "p^ = 0" approximation 
simplifies to 

= approximation, 3x3 mixing: Zi = — , (23) 

1 + 4 + 

as a direct generalization of eq. (22). 

The matrix Zn defined through eqs. (20)-(21) goes over into a unitary matrix Rn in this 
approximation, 

(Rll i?i2 Rl3\ 
R21 R22 R23 I • (24) 
^31 R32 R33 / 

The matrix Rn diagonalizes the matrix Mn(0) arising from eq. (4) in the "p^ = 0" ap- 
proximation. Rn can therefore be used to connect the mass eigenstates /ii, /12, with the 
original states h, H, A, 

'h^\ /h\ (Mlr>^=o 

hA = Rn • U , RnMn(0)Ri = M^^^.^^ « | . (25) 
W,.=o Vl V Ml^,._ 
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Since Mn is hermitian, the matrix Rn is unitary. 

As shown in Ref. [16], a better approximation of the full result can be achieved by 
defining the effective couplings in the "p^ on-shell" approximation. The unitary matrix Un 
is defined such that it diagonalizes the matrix Rc (Mn(]5^ on-shell)) arising from eq. (4) 
in the "p^ on-shell" approximation and neglecting the imaginary parts. This yields 

p2 on-shell approx.: I /i2 I = U„ ■ I I , Un = I C/21 i722 C/23 I , (26) 



p2 on— shell 




on-shell 

Un Re (Mn(p' on-shell)) Ut = M^^^, 

V 2 , ,, , 

\ 'i3)P on— shell/ 

The elements of Un can be used to quantify the extent of CT'-violation. For example, 
C/fg can be understood as the CP-odd part in hi, while + UI2 make up the CP-even 
part. The unitarity of Un ensures that both parts add up to 1 . 

The elements of Un can be interpreted as effective couplings of internal Higgs bosons in 
a loop diagram, which take into account leading higher-order contributions. As an example 
we show the couplings of a neutral Higgs boson to two gauge bosons, VV = ZZ, W^W~ . 
Beyond the lowest order in the cMSSM all three neutral Higgs bosons have a CP-even 
component, so that aU three Higgs bosons have non-vanishing couplings to two gauge 
bosons. The couplings normalized to the SM values are given by 

9hiVV = Uii sin(/3 - a) -|- Ui2 cos(/3 - a) . (27) 



4. Conclusions 



Large higher-order corrections enter the Higgs boson sector of the MSSM via Higgs-boson 
self-energies. Their effects have to be taken into account for the correct treatment of loop- 
corrected Higgs-boson mass eigenstates as external (on-shell) or internal particles in loop 
diagrams. We have shown how the loop corrections can be taken into account in the treat- 
ment of external (on-shell) Higgs bosons, and how effective couplings, taking into account 
leading higher-order effects, can be constructed. 

The proceedure presented here is implemented in the code FeynHiggs2.5 
[16,27-29]. For the Higgs-boson self-energies the complete one-loop corrections 
[16,30-32] are supplemented by the available two-loop corrections in the Feynman- 
diagrammatic approach for the MSSM with real parameters [27,28,33] and a resummation 
of the leading (s)bottom corrections for complex parameters [34]. Higgs boson decays 
and production cross sections [35] are evaluated using Zn (20,21), ensuring the on-shell 
properties of the external particles. Effective Higgs boson couplings are obtained with the 
help of Un (27). In this way FeynHiggs provides a precise prediction of Higgs boson 
production and decay properties and the respective coupUngs. 
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